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We present direct representations of the scaling functions of the 3d 0(A) model which are rele- 
vant for comparisons to other models, in particular QCD. This is done in terms of expansions in 
the scaling var iablez = f"//i 1 /' 35 . The expansions around z = and the corresponding asymptotic 
ones for z — > ±°° overlap such that no interpolation is needed. We explicitly present the expan- 
sion coefficients which have been determined numerically from data of a previous high statistics 
simulation of the 0(A) model on a three-dimensional lattice of linear extension L = 120. This 
allows to derive smooth representations of the first three derivatives of the scaling function of the 
free energy density, which determine universal properties of up to sixth order cumulants of net 
charge fluctuations in QCD. 
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1. Introduction 

We provide representations of the scaling functions of the three-dimensional 0(4) model 
which can be used in tests of other models on their membership of the corresponding universality 
class. Contrary to the often analyzed scaling function of the order parameter, the so-called mag- 
netic equation of state, our main interest here is to determine directly the scaling function of the 
free energy density and its derivatives. This is especially of importance for applications to quantum 
chromodynamics (QCD) with two degenerate light-quark flavors at finite temperature. Two-flavor 
QCD is believed [jl|]-[§] to belong to the 3d 0(4) universality class at its chiral transition in the 
continuum limit. In the vicinity of the chiral phase transition temperature the reduced temperature 
variable in QCD also depends quadratically on the quark chemical potential. Derivatives of the 
singular part of the free energy density of QCD with respect to chemical potential, which define 
cumulants of fluctuations of net quark number, thus are controlled by scaling functions that are 
given by derivatives of the scaling function of the free energy density in a three-dimensional 0(4) 
model. 

Obtaining explicit parametrizations of higher order derivatives of the scaling functions of the 
free energy density became of interest recently as these higher order derivatives control the scaling 
behavior of fluctuations of conserved charges, e.g. the net baryon number [Q]. These quantities are 
currently measured at RHIC []|] and will also be measured in heavy ion experiments at the LHC. 

2. The three dimensional 0(4) model 

The specific model which we study here is the standard 0(4) -invariant nonlinear a-model, 
which is defined by 

<x,y> x 

where x and y are nearest-neighbor sites on a three-dimensional hypercubic lattice, and fa is a 
four-component unit vector at site x. The coupling J and the external magnetic field H are reduced 
quantities, that is they contain already a factor /3 = l/T. In fact, we consider in the following the 
coupling directly as the inverse temperature, / = l/T. The partition function is then 

Z(T,H) = Jl[d%8($i-l)ex V (-pj?) . (2.2) 

X 

We introduce the order parameter M as the derivative of the free energy density, f(T,H) = -^lnZ, 
with respect to the magnitude of the external magnetic field H = Hen, 

M=-^ = (^), (2.3) 

where $H is the field component parallel to the magnetic field H. 

In the vicinity of the critical point the free energy density may be splitted into a singular (non- 
analytic) (f s ) and a non-singular (J ns ) part, 

f(T,H) = f s (T,H)+f ns (T) . (2.4) 
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The singular part is a homogeneous function of the variable h = H/Ho and the reduced temperature 
t = (r — T c )/Tq, where Hq and 7b set the scale in the critical region. The singular part may be 
expressed in terms of a universal scaling function //, which itself only depends on the scaling 
variable z = t/h l l A , i.e., 

f s = H h l+ ^ s f f (z) . (2.5) 

Here we have introduced the gap exponent, A = which is given in terms of the more commonly 
used critical exponents j8 and 8. The latter define the scaling properties of the order parameter as 



function of temperature at h = and as function of the external field at t = 0, respectively. Eq. 2.5 



establishes the relation between the universal scaling function of the order parameter (fa) and the 



scaling function of the free energy density (//). Using Eq. ^3J we find 

M = h l ' s f G {z) (2.6) 

r 



fG{z) = -(l + \)f f {z)+ Z -f' f {z). (2.7) 



In the following we will exploit the differential equation, Eq. |2.7[ , to determine the scaling function 
f f (z) from f G (z). 



3. Scaling functions of the free energy density and the order parameter 



As the universal scaling functions of the free energy density, ff(z), and the order parameter, 
fc{z), are related through the differential equation, Eq. ^7^, the knowledge of fc(z) is sufficient to 
determine //(z). We summarize in the following the relevant relations that determine f/(z), once a 
suitable parametrization of fo{z) is known. Further details are given in Ref. [^]. 

We consider a parametrization of fc{z) by introducing three series expansions that are valid 
for small z and in the asymptotic regions z — > ±°°, respectively, 



fciz) 



^2n=0 bnZ U 



7-7 



-2nA 



^2n=Q^nZ 



, z small 

, z +°° 
, z -> — °° 



(3.1) 



The corresponding parametrization for the scaling function of the free energy density is then 
given by, 

YZ=o a nZ n , z small 



//(*) 



,2- a 



J2n=0 c n z 



-2nA 



-)-oo 



(3.2) 



The relation between the expansion coefficients in the series representations for fa and fj are easily 



obtained by using the differential equation, Eq. |2.7| , and comparing coefficients for n > 0, 

Ab„ , -d+ _ 2d~ 



a + n — 2 



2(n + l) 



n + 2 



(3.3) 
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bo 


bx 


b 2 


b, 


1 


-0.3166125 ±0.000534 


-0.04112553 ±0.001290 


0.00384019 ±0.000667 








K 




0.006705475 ±0.001704 


0.0047342 ±0.001429 


-0.001931267 ±0.000312 




K 








0.007100450 ±0.000160 


0.0023729 ±0.000095 


0.000272312 ±0.000021 



Table 1: Coefficients of the small z-expansion of the scaling function /g(z) of the order parameter. For 
n > 4 we give different expansion coefficients for negative and positive z-values. 

and Cj~ = 0. This leaves the coefficients Cq still undetermined. They can be obtained as 

A f'°° 

c o = o-a / dyy "~ 2 ifo(y)-fo(o) -yfM] , 0.4) 

* — t* Jo 

c o = f d y (-y^ 2 ifa(y)-fa(P) -y/g(0)] • (3-5) 

Here a is the specific heat critical exponent, which is negative in the three dimensional 0(4) 
universality class. 

The coefficients a n for all n and c„, for n > 0, are obtained from a parametrization of the scaling 
function of the order parameter. The relevant expansion coefficients b„, b^ and d^ have been 
obtained from numerical data for the order parameter itself as well as its susceptibility [Q]. In this 
step one explicitly makes use of a set of values for critical exponents in the three dimensional 0(4) 
universality class. We used: j8 = 0.380 and 8 = 4.824. All other critical exponents can be derived 
using hyperscaling relations. E.g., the specific heat exponent is given by, a = —0.213. We list the 
resulting expansion coefficients in Table [l] and ^. Note that we give different expansion coefficients 
b„ and b~ for n > 4 to better reproduce the asymmetric form of the scaling function fc(z) also for 
small values of z with only a small number of expansion coefficients. The corresponding scaling 
function of the order parameter and its first derivative is shown in Fig. |]. 

Having at hand a parametrization of fc{z) we finally can determine the remaining coefficients 
Cq , which complete the parametrization of //(z). These expansion coefficients are listed in Table ||. 



4 


4 


d+ 


1.10599 ±0.00555 


-1.31829±0.1087 


1.5884 ±0.4646 


d 


d- x 


d 2 


1 


0.273651 ±0.002933 


0.0036058 ±0.004875 



Table 2: Coefficients for the asymptotic series expansions of /g(z) in the region of large positive and 
negative z-values, respectively. 

4. Discussion and Conclusions 

The availability of high accuracy numerical data on the order parameter and its susceptibility 
in a three dimensional, 0(4) symmetric spin model allowed us to extract the underlying scaling 
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Figure 1: The scaling function of the order parameter (left) and its derivative (right). 
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Figure 2: The scaling function of the free energy density and its first three derivatives. 
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4 


c o 


0.422059886 ±0.010595 


0.229176194 ±0.010669 



Table 3: The leading expansion coefficients for the singular part of the free energy density. 

function of the free energy density and its first three derivatives. As the specific heat exponent a is 
negative in the 3d 0(4) universality class, it is only the third derivative with respect to temperature, 
which diverges at the critical point. The corresponding scaling function f" (z) has two extrema; a 
rather shallow minimum in the symmetry broken phase and a pronounced maximum in the sym- 

3 

metric phase. The latter is located at Zp — 1.45. This happens to be close to the location of the 

2 

peak in the susceptibility of the order parameter, Zp = 1.374(3). 

The higher order derivatives of the scaling function of the free energy density play a central 
role in the discussion of fluctuations of conserved charges in QCD, e.g. the singular behavior of 
the 2«-th order cumulant of net baryon number fluctuations is related to the n-th derivative of f/(z). 
The change of sign of ff'{z) and its pronounced maximum characterize the QCD transition. In fact, 
the change of sign of ff(z) suggests that 6th order cumulants of net baryon number are negative 
in the vicinity of the QCD transition line. This may be detectable in a heavy ion collision, if the 
production of hadrons (freeze-out) occurs at temperatures and baryon chemical potentials that are 
close to the QCD crossover transition line. 
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